Charge ordering in the spinels AIV2O4 and LiV204 
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We develop a microscopic theory for the charge ordering (CO) transitions in the spinels AIV2O4 
and LiV2 04 (under pressure). The high degeneracy of CO states is lifted by a coupling to the 
rhombohedral lattice deformations which favors transition to a CO state with inequivalent V(l) 
and V(2) sites forming Kagome and trigonal planes respectively. We construct an extended Hub- 
bard type model including a deformation potential which is treated in unrestricted Hartree Fock 
approximation and describes correctly the observed first-order CO transition. We also discuss the 
influence of associated orbital order. Furthermore we suggest that due to different band fillings 
AIV2O4 should remain metallic while LiV2 04 under pressure should become a semiconductor when 
charge disproportionation sets in. 

PACS: 71.30.-|-h; 71.70.-d; 71.10.Fd; 71.20.Be 

It was Wigner [1] who first pointed out that electrons 
form a lattice when the mutual Coulomb repulsion dom- 
inates the kinetic energy gain from delocalization. He 
considered an electron gas with a positive uniform back- 
ground (jellium). Within that model a lattice will form 
only when the electron concentration is very low. De- 
tailed calculations have determined the critical concen- 
tration below which lattice formation or charge ordering 
takes place and a value of ~ 35aB was found [2] . Here 
Tc is the critical value of the average distance between 
electrons which must be exceeded and as is Bohr's ra- 
dius. The conditions for charge ordering of electrons are 
much better when the uniform, positive background is 
replaced by a lattice of positive ions. Depending on the 
overlap of atomic wavefunctions on neighbouring sites 
the energy gain due to electron delocalization can be 
rather small and a dominance of Coulomb repulsion is 
more likely. For a more detailed description see, e.g. [3]. 
A good example is Yb4As3 which exhibits charge order 
(CO) and for which a microscopic theory was provided 
[4]. Here Yb 4/ holes gain a small energy only by de- 
localizing via As 4p hybridization. Another well known 
example is magnetite Fe3 04, a spinel structure. Verwey 
and Haayman [5] found a metal-insulator phase transi- 
tion to occur which they attributed to a charge ordering 
of the Fe^+ and Fe'^+ sites on the pyrochlore sublattice of 
the spinel lattice (B-sites) . According to Verwey [5] in the 
charge ordered state the corner-sharing tetrahedra form- 
ing the pyrochlore structure are occupied by two Fe^+ 
and two Fe"^^ ions each. This rule, (tetrahedron rule) was 
most clearly formulated by Anderson [6] and allows for an 
exponential number of different configurations. Recently 
this rule has been questioned by LDA+U calculation [7] 
which seemingly give good results for the observed charge 
disproportionations [8] when the low-temperature struc- 
ture is put into the calculations. 

In this investigation we consider AIV2O4 [9] and 
LiV204 [10,11]. In both systems the V ions form a py- 



rochlore lattice. While in AIV2O4 a charge order phase 
transition at approximately 700 K has been observed 
[9] at ambient pressure, in LiV204 electronic charge or- 
ders only under hydrostatic pressure of approximately 6 
GPa [11]. In both cases structural changes are associ- 
ated with charge ordering. In AIV2O4 the lattice dis- 
tortions change the cubic high-temperature phase into a 
lattice with alternating layers of Kagome and triangular 
planes. They are stacked in [111] direction and formed 
by nonequivalent V(l) and V(2) sites. The rhombohe- 
dral lattice deformation causes an increase in the V(l)- 
0(1) and V(l)-0(2) bond lengths and a decrease in the 
V(2)-0(l) bond length. The changes in the distances 
V(l)-V(l) denoted by Li and V(l)-V(2) by L2 below the 
phase transition temperature Tco define a deformation 
e(T) = (L2 — Li)/L. L is the V-V distance in the undis- 
torted case. At T = 300 K this value was found to be 
e{300K) = 0.016. Below Tco the average valence of the 
V ions changes from -t-2.5 in the high-temperature phase 
to -\-2.5-5 at the V(l) sites and to -f 2.54-3(5 at the V(2) 
sites. This is schematically shown in Fig. 1. 
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FIG. 1. (a)Sublattice of the V ions. In the presence of a 
distortion the V(l) sites form a Kagome and the V(2) sites 
a triangular lattice, (b) AIV2O4 and LiV204 in the atomic 
limit. Shown are the splittings of the t2g orbitals, their oc- 
cupations, the spin vectors and the angle a with respect to 
the [111] plane. 

Also associated with the charge-order transition are 
a small but steep increase in the resistivity and a pro- 
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nounccd decrease of the magnetic susceptibility. In 
LiV204 the lattice distortion connected with the charge 
order under pressure seems to be similar to that of 
AIV2O4 as powder X-ray diffraction indicate [11]. How- 
ever the structural parameters have not yet been deter- 
mined. Due to the different filling of the t2g band AIV2O4 
seems to remain a semimctal or a small gap semiconduc- 
tor while LiV204 is becoming an insulator with a steep 
increase in the resistivity below the phase transition tem- 
perature Tco- 

The first thought is to perform for AIV2O4 and LiV204 
under pressure LDA+U calculation of the type reported 
for FC3O4 (see [7]). We have done such calculations for 
AIV2O4. However, although the LDA-I-U approach was 
shown to provide a realistic description of the electronic 
structure of magnetic insulators, its application to para- 
magnets is less justified because calculations can only 
be performed for a magnetically ordered state. In addi- 
tion, LDA-I-U is too crude an approximation to reproduce 
the low energy excitation spectrum of strongly correlated 
metals and, thus, can hardly be used to study the high 
temperature metallic phase of AIV2O4 and the change of 
the electronic structure upon heating through the phase 
transition. For LiV204 LDA-I-U calculations are not yet 
possible since the CO crystal structure under pressure is 
not known. Therefore we proceed here differently. We 
want to provide a microscopic Hamiltonian in order to 
identify the physical process which results in the observed 
structural phase transition and the accompanying charge 
order while enforcing a proper nonmagnetic state. The 
price are adjustable parameters entering the theory for 
which no ab initio values are available. From standard 
band-structure calculations it is known that the Fermi en- 
ergy Ep of V 3d electrons lies within the t2g band which 
is well separated from the remaining bands. Therefore 
we only include t2g electrons in the model defined by 

H = Ho + Hint + He-p , with (1) 
Hint = '^{{U + 'iJ)'^nii^^ini^^i + U'^ ni^^^ani^^'a 
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Here Hq describes the kinetic energy. The electron cre- 
ation and annihilation operators are specified by four in- 
dices, i.e., for the unit cell (denoted by /), the sublat- 
tice (/X = 1,2,3,4), the t2g orbital = d^y, dy^, dzx)i 
and the spin {a =T)i)- The brackets (...) and ((...)) in- 
dicate a summation over nearest-neighbour (n.n.) and 
next-nearest-neighbour (n.n.n.) sites respectively. The 



term Hint describes the on-site Coulomb and exchange 
interactions U and J among the t2g electrons. The last 
term contains the Coulomb repulsion V of an electron 
with those on the six neighbouring sites. Finally Hep 
describes the coupling to lattice distortions. The defor- 
mation potential is denoted by A, it is due to a shift in 
the orbital energies of the V sites caused by changes in 
the oxygen positions relative to the vanadium positions. 
While the energy shift is positive for the V(l) sites it is 
negative for the V(2) sites. 

The elastic constant K refers to the C44 mode and de- 
scribes the energy due to the rhombohedral lattice de- 
formation. It is reasonable to assume that like in Fe304 
[12] and Yb4As3 [13] only the C44 mode is strongly cou- 
pled with the charge disproportionation. We can give 
at least approximate values for all parameters except V 
and A. Their ratio will be fixed by the CO transition 
temperature while keeping the constraint V<CU. For the 
on-site Coulomb- and exchange integrals we take U = 3.0 
eV and J = 1.0 eV which are values commonly used for 
vanadium oxides [14]. Band structure calculations which 
we have performed demonstrate that the hopping ma- 
trix elements t'','',, between different orbitals v ^ v' are 
negligible. For simplicity we first ignore them and omit 
v^v' . Then t = —t^^j_i{l,l') when U I' are n.n. and t' for 
n.n.n.. We will take into account the orbital dependent 
hopping matrix elements coming from a tight-binding fit 
of LDA calculations as orbital order is included. Further- 
more the C44 elastic constant is not known for AIV2O4 
or LiV204, while computational methods for its ab initio 
calculation in the case of materials with strong electronic 
correlatations are not mature enough. Therefore we use 
a representative value 044^^^ = 6.1 • 10^^ erg/cm^ for 
AIV2O4 where Vl is the volume of the cubic unit cell with 
a lattice constant of a = 5.844 A. This value is close 
to the experimental value for Fe304 which has also the 
spinel structure. This leads to ~ 1.1 ■ lO^eU. The 
deformation potential A is not known but is commonly 
of the order of the band width. For convenience we in- 
troduce the dimcnsionless coupling constant A = ^ jKt 
and lattice distortion 5^ = eA/t. From LDA calculations 
the bandwidth is 8t=2.7 eV, and therefore a reasonable 
value is At = A^/JsT = IcV. This means A = lO.SeV 
which is twice the value of Yb4As3 [4]. We begin by con- 
sidering Hq only. By transforming to Bloch states aX^^„ 
we obtain 



£1,2 (k) = 2{t- t') 



(2) 



,4{k) = -2t 1 T + % +2i' l-2ryk±2v^ 



Here ryk = (cos fca, cos fcy 



cos ky cos kz 



■ cos kz cos kx)- 



With four V ions per unit cell, three t2g orbitals and two 
spin directions there are altogether 24 bands, of which 
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twelve are dispersionless and degenerate (€1,2)- In ad- 
dition there are two sixfold degenerate dispersive bands 
present (£3,4). The band structure has been considered 
before [15] and we refer to that reference for more details. 
The Fermi energy lies in a region of high density of states 
(DOS) N(E) which depends sensitively on the ratio t' /t. 

Next we discuss the effect of the deformation-potential 
coupling described by Hep- When it is sufficiently strong 
it may lead to charge order. However, there is no gap 
opening but only a sharp decrease of DOS around Ej?, 
i.e., the system remains metallic. When we include the 
U and V terms in mean-field approximation, assuming a 
paramagnetic state, the results remain qualitatively un- 
changed. A strong on-site U term suppresses CO while 
increasing inter-site V term will induce CO. Again there 
is no gap opening at Ej? but a decrease of DOS at the CO 
transition for all t' . We note that in the present case a 
n.n.n hopping t' supports CO because it increases N(Ef), 
while it prevents CO in the half-filled square lattice [16]. 
Since here t' « t has little effect on CO we neglect it 
for simplicity. 

In our calculations we assume a homogeneous defor- 
mation along [111] direction. This yields the following 
relation between the lattice distortion and charge dispro- 
portionation 



A (^2 - ni) /2 



(3) 



where ni = X)iy(T('^iii'<^) is the occupation of a V(2) site 
while 71-2 = ^3 = n4 with = ^^„{ni^vrj) is the occu- 
pation of a V(l) site. 

In any case, the mean-field analysis of the model 
Hamiltonian (1) leads to a second order phase transi- 
tion as function of temperature. This is at odds with 
the experimental observation of a strong first order CO 
transition. We believe that this is the effect of charge cor- 
relations caused by large U and V. To include this effect 
in a qualitative way, we have treated the model within 
an unrestricted Hartree-Fock calculation by breaking the 
magnetic symmetry but preserving the constraint of zero 
total moment. We assume that the sites ^ have an oc- 
cupation = '^^crinifiiya) and a magnetization = 
'^^^o'ini^va) ■ The spins are assumed to be directed to- 
wards the center of the tetrahedron in the undistorted 
phase. In the distorted phase we change the angle a so 
that the net magnetization remains zero (see Fig. 1 (b)). 
Because there are two tetrahedra per unit cell transla- 
tional symmetry is maintained. The free energy is 
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where N is the total number of unit cell and A^e^lO N is 
the electron number per unit cell. Furthermore S is the 
grand canonical partition function 



InS = 3^ In 



-/3(e^(k)-Jl) 



(5) 



with (3 — l/ksT. The chemical potential fi is adjusted 
to yield the correct filling, i.e., Ng — ^-^InE. The en- 
ergy bands e|^(k) depend now on spin a for which the 
[111] direction is chosen as quantization axis. They are 
computed by diagonalizing the Hamiltonian (1) in mean- 
field approximation and by self-consistent determination 
of and from (9F/9np)— = and [dF / dm^)— = 0, 
respectively at various temperatures T. 

With the magnetization pattern described above and 
shown in Fig. 1 (b) we must transform the nearest- 
neighbour hopping matrix elements so that the different 
spin directions arc accounted for. We label the electronic 
spin functions of sublattice 1 by \ai) and |cti) and those 
referring to sublattice 2, 3, and 4 by \(jf_i) and We 
then apply the transformation to a local spin quantiza- 
tion axis 
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(6) 



where 6* = f -h a and 7^, = 0, ^, ^ for ^i=2, 3, 4. All 
the interaction terms as well as the electron-phonon cou- 
pling term remain invariant under this transformation. 
The angle a is adjusted such that m^—Q for all dis- 
tortions. 




FIG. 2. Charge disproportionation as function of V/t. 
The rifi denote the occupation numbers of the four sites of a 
tetrahedron. In the inset the changes in the density of states 
are shown when V/t is just below and above the critical value 
at which charge ordering sets in. At the critical V/t=1.67 
m =2.5-1-3(5 and ni=2.5-5 (i=2-4) with a charge dispropor- 
tionation S ~ 0.25. Here U=3.0eV, J=1.0eV, \t=l.GeV and 
8t=2.7eV. 

The result for the charge disproportionation at T—Q 
as function of the electronic nearest-neighbour Coulomb 
repulsion V is shown in Fig. 2. As usual the charge dis- 
proportionation is larger (by a factor of 2.5 here) than 
the one obtained from a valence bond analysis of the dis- 
torted structure. This is due to the simplified Hamilto- 
nian which cannot describe the screening effects of non-d 
electrons. The same observation was made for Yb3As4 
[4]. Also shown as an inset is the change in the DOS for a 
V value just below and above the critical value at which 
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charge order does occur. Apparently the DOS changes 
considerably near Ep due to charge order. This may 
explain the drop in the susceptibility and the small but 
steep increase in the resistivity in AIV2O4 when charge 
order sets in. The self-consistent field calculations were 
done for different temperatures. The phase transition to 
the distorted phase is found to be of first order. Fix- 
ing the value V/t=1.67 leads to the calculated transition 
temperature is Tc=660 K in reasonable agreement with 
the experimental value of T®'^p=700 K. The phase dia- 
gram is found to be very simple. For U=3 eV, J=1.0 eV 
and 8t=2.7 eV we find a phase boundary line of the form 
Xt+SV=7.8t. For all Xt values above that line the trigo- 
nal, i.e., distorted phase is found while for values below 
that line the phase remains cubic. 

It is obvious that with a lattice distortion an associ- 
ated orbital order will occur [17]. The point symmetry 
in the charge ordered state is Dsa and lifts partially the 
3-fold degeneracy of the t2g states. This results in one 
aig orbital and 2-fold degenerate orbitals. Experi- 
ments show in the distorted phase an elongation in [111] 
direction and a corresponding contraction perpendicular 
to that direction. Therefore the orbital energy of the aig 
state is expected to be higher than that of the e'g states. 

From the behavior in the atomic limit (see Fig. 1 (b)) 
we expect that the system remains gapless in the charge 
ordered state even when orbital order occurs. This is in- 
deed what we find. We have adopted hopping matrix el- 
ements which reproduce the LDA band structure. They 



opening of a gap is in agreement with the behavior of 
p{T) as found in Ref. [11]. 



are to = t^^f '"'"=-0.525 eV and ii 
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0.152 eV. The upmost band remains almost fiat. By solv- 
ing self-consistently for the mean-field parameters n^i/ 
and nifii, under the constraint that the occupation num- 
bers of the three V(l) sites remain equal we find that 



— TI2" ^11^2 1 "'S — "'3 ^ ""3 cLiita /fc^ — tt^ ^11^ 

For simplicity we have used here the notation of the t2g 
basis. It shows that the orbital orders on sublattices 2, 3 
and 4 are perpendicular to each other. At the V(2) sites 
there is no orbital ordering in the charge ordered phase. 
Since the system remains gapless and the orbital order 
does not change the effects induced by charge order and 
magnetic order, i.e., a sharp decrease of DOS around the 
Fermi surface and a first-order phase transition respec- 
tively, its role is insignificant at 5/12-filling (AIV2O4). 

Now we compare the above findings with the charge 
order observed in LiV2 04 under pressure. Here the av- 
erage d electron number per V site is 1.5. The atomic 
limit shown in Fig. 1 (b) suggests a gap and therefore 
insulating behavior in the charge ordered phase. This is 
also what is found when the Hamiltonian (1) is treated 
in mean-field approximation for 1/4 filling as required 
for LiV204 instead of 5/12. The essential role for or- 
bital order and gap formation is played by the different 
3d-occupation of LiV204 therefore we use similar values 
for hopping terms and interaction strengths are used as 
for AIV2O4. The DOS at T=0 is shown in Fig. 3. The 
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FIG. 3. Density of states or *'LiV2 04 in tlic distorted 
phase. It is soon that the system is a semiconductor 
with a gap £;5=0.49oV. Here U=3.0eV, J=1.0eV, V=0.7eV, 
AVA'=1.0cV, i(,=-0.53cV, fi=0.15eV 

In summary, we have provided a microscopic model for 
the observed charge order in AIV2O4 and LiV204 under 
pressure which is mainly driven by a deformation po- 
tential coupling to t2g states. The model accounts cor- 
rectly for the observed first order CO transition. We 
have also shown that due to a difference in the filling 
factors AIV2O4 should remain metallic in the charge or- 
dered phase while LiV204 should become an insulator or 
semiconductor. 

We thank Drs. G. Venketeswara Pai, T. Maitra, E. 
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